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Abstract 

We have developed two analytical models to describe the performance of cryo- 
genic microcalorimeters and bolometers. One of the models is suitable to describe 
Transition Edge Sensor(TES) detectors with an integrated absorber, the other is 
suitable for detectors with large absorbers. Both models take into account hot- 
electron decoupling and absorber decoupling. The differential equations describing 
these models have been solved using block diagram algebra. Each model has pro- 
duced closed form solutions for the detector's responsivity, dynamic impedance, and 
noise equivalent power for phonon noise, Johnson noise, amplifier noise, 1/f noise, 
and load resistor noise. 
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1 Introduction 



The operational principle of simple microcalorimeters and bolometers is based 
on three components. An absorber where the incoming power or energy is 
dissipated and converted into a change in temperature, a sensor that reads 
the change in temperature, and a thermal link from the detector to the heat 
sink that brings the system back to equilibrium after a measurement. The 
sensor is usually a resistor whose resistance depends strongly on temperature 
around the working point. In this change in resistance can be measured 

as a change in voltage or a change in current using a current or voltage bias. 

In 1982 J.C Mather [1] presented a complete non-equilibrium theory for the 
noise in simple bolometers with ideal resistive thermometers and in 1984 it was 
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extended to microcalorimeter performance [2]. At temperatures below 200 mK 
the ideal assumptions are no longer valid and complex thermal architectures 
are needed to understand the behavior of these devices. At these low tem- 
peratures the thermal fluctuations between the thermometer lattice and its 
electron system, or between the thermometer and the absorber are no longer 
negligible and therefore these components must be considered as separate en- 
tities in the model. These non-ideal effects are called electron decoupling and 
absorber decoupling. Another consequence of working at low temperatures is 
the increased dependence of the thermometer resistance on the readout power, 
making the ideal resistance-temperature relationship inaccurate. Galeazzi and 
McCammon [3] constructed a general procedure for developing bolometer and 
microcalorimeter models for these complex thermal architectures using block 
diagram formalism of control theory. 

To quantify the relation between incoming power or energy and the measured 
change in voltage or current, including non-ideal effects, this paper follows 
the modeling procedure of Galeazzi and McCammon [3] . The first step in this 
modeling procedure is to set up the temperature equations, then apply a Tay- 
lor expansion to derive a linear model for small temperature deviations from 
equilibrium. Afterward, Fourier transforms are used to express the equations 
in the frequency domain, and finally the coupled equations are solved using 
block diagram algebra. This procedure yields closed form solutions for the 
responsivity and dynamic impedance of the model, including noise contribu- 
tions. 

The two models developed in this paper are modifications of the absorber de- 
coupling model obtained in [3]. Model 1 was developed to describe the new 
generation of transition edge sensor detectors where the absorber is not elec- 
trically isolated from the thermometer by a gluing agent but rather the two arc 
deposited one on top of the other [4]. Model 2 describes microcalorimeters that 
have the heat sink connected to the absorber instead of to the thermometer. 
This may occur when the absorber is much bigger than the thermometer and 
therefore it is necessary to connect the heat sink to the absorber rather than 
to the thermometer. These two models will help optimize the next generation 
of detectors, and because of the analytical results of the modeling procedure 
the relations between the detector's resolution and the different parameters 
included in the model should be clear. 



2 Model 1 



This model is suitable to describe TES detectors with an integrated absorber. 
In this model both the electron system and the absorber are detached from the 
lattice. The lattice is connected to the heat sink by a thermal conductivity G, 
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the electron system is connected to the lattice by a thermal conductivity G e -i, 
and the absorber is connected to the electron system by a thermal conductivity 
G a ( see Fig. 1). The absorber is directly connected to the electron system 
instead of the lattice because with integrated absorbers the absorber-lattice 
thermal coupling is expected to be negligible compared to that of the absorber- 
electron system. 



2.1 Responsivity S(uj) 



The following equations determine the temperature for each of the three com- 
ponents in the model: 



C a ^ + jG a (T>)dT' = W (1) 



d{T' e ) 



rp/ rp/ 



Ce ^t + / GaiT '^ dT ' + / G e-l( T ') dT> = P ( T e) ( 2 ) 



T; 

d(T{) 



Cl dt + I G *- l( - T "> dT ' + / G(T')dT' = 0, (3) 



where C a , C e , and C% are the heat capacities of the absorber, the electron 
system, and the lattice system respectively, and T' a , T' e , and T( are the corre- 
sponding temperatures. W is the incoming outside power to be measured and 
P{T' e ) is the Joule power dissipated into the sensor by the bias current /volt age. 
In the case of microcalorimeters W = ES(t ), where E is the photon energy 
and 8 {t ) is the delta function. 

The equilibrium conditions of the system are obtained by setting the outside 
power to zero {W = 0), and d{T' x )/dt = (x— a, e, or I) since the equilibrium 
temperatures are independent of time. Therefore the equilibrium temperatures 
T a of the absorber, T e of the electron system, and TJ of the lattice are given by 
the integrals in the previous three equations. For example the integral in Eq. 1 
must equal zero at equilibrium, which implies that the thermal equilibrium 
temperature of the absorber is the same as that of the electron system. We are 
interested in small deviations about the equilibrium temperatures, therefore 
we set T' x = T x + AT X , where T x is the equilibrium temperature for each 
component of the model, and AT X is the small temperature deviation from 
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equilibrium: 

Ta+ATt 



Ca dt + l G ^ dT ' = W W 

T e +AT e 

d(T + AT ) Te+ r ATe Te+ ? Te 
Ce { e dt e) + J G a (T')dT'+ j G e ^(T')dT> = P(T e + AT e )(h) 



Ci diTl+ d ^ Tl) + / G e _ l (T')dT'+ J G(T')dT' = 0. (6) 

T e +AT e 



Ta+ATa Ti+ATi 
T ( +AT ; T ; +AT ; 



In the small signal limit AT X is small compared to T x , and a Taylor expansion 
up to the first AT X term is appropriate. The results are the equations that 
determine small temperature deviations about equilibrium: 

Ca ^^~ + ° aATa = W + ° aATe (7) 
Ce^ft + G a AT e + G e _ t (T e )AT e = AP + G a AT a + G e _,(T,)A7i (8) 

+ G e -i(Ti)ATi + GAT; = G e -z(T e )AT e , (9) 



where AP = P(T e + AT e ) - P(T e ), and for simplicity we used G a = G a (T a ) 
and G = 

These are coupled differential equations which are difficult to solve directly; 
instead they are transformed into coupled algebraic equations using Fourier 
transforms. The quantity AP represents what is known as the electro-thermal 
feedback term and can be written as AP = — GetfAT b , where Getf — 
P(R — Ri)a/T e R(RL + R); (see reference [3]). Converting Eqs. 7, 8, and 9 
into the frequency domain using Fourier transforms we obtain: 

jujC a AT a + G a AT a = W + G a AT e (10) 

juC e AT e + (G a + G e -i(T e ) + G ETF )AT e = G a AT a + G e -i(Ti)AT t (11) 

juCiATi + (G e _,(T,) + G)AT t = G e ^(T e )AT e . (12) 



With these equations it is possible to solve for AT e , which can be related to 
the measured quantities AI or AV, using the typical detector readout circuit 
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of Fig. 2: 



AV=v iwh AT ' (13) 
A/ -- 7 tra Ar " (14) 



where a = T e /R x dR/dT e is the sensitivity of the detector, R L is the load 
resistance, R is the resistance of the detector, V is the voltage across R, and 
/ is the current flowing through R. 

To simplify the notation let X be either V or /, and introduce the quantity 
A tr = R/X x dX/dR to be deduced from the previous two equations. Then 
Eqs. 13 and 14 can be summarized as: 

— = «Ar^- (15) 



Equations 10, 11, 12, and 15 can be solved using the block diagram of Fig. 3. 
To set up the block diagram consider the left hand side of Eqs. 10, 11, 12, and 
15 as the response function of the absorber system, electron system, lattice 
system, and circuit readout respectively. The right hand side of these equations 
corresponds to the input to each system. Connecting the response functions 
with their appropriate inputs leads to the block diagram of Fig. 3. 



To solve the block diagram in Fig. 3 for Al(w) we used the procedure and 
simplification rules of the block diagram formalism described in [3] . This result 
is then used to find the responsivity, which is defined as S(u) = AX(u)/W(u). 
The following responsivity characterizes the response of Model 1 detectors: 



S(w) = 



(G a + JLuC a ) 



(G a + G e -i(T e ) + Getf + juC e ) 



G e -i(Ti)G e -i(T e ) 



XaA 



-Gl 



^.(16) 



2.2 Dynamic Impedance 



A detector can also be described by its complex dynamic impedance Z(uj) = 
dV{uj)/dI{uj). The dynamic impedance differs from the detector resistance due 
to the effect of the electro-thermal feedback. When the current changes, the 
power dissipated into the detector changes too, therefore the temperature and 
the detector's resistance change. The dynamic impedance is a useful parameter 
because it is easily measured experimentally. To find the dynamic impedance 
we use -G ETF AT e = AP in Eqs. 11, and use Eqs. 10, 11, and 12 to find AT e 
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in terms of AP, u, the heat capacity of each of the three components, and the 
three thermal conductivities: 



dP 

dT e = - , . (17) 

[juC e + G a + G e -t(T e ) - G J^ Ca - fjc+aZlT^+G / 



Differentiating Ohm's law (V = IR) and using the definition of sensitivity a 
we obtain: 

dV = RdI + I^^. (18) 



Substituting Eq. 17 into Eq 18 and using the fact that dP = Vdl + IdV, it 
is possible to solve for dV/dl and obtain the following result for the dynamic 
impedance: 



°P + T[juC e + G a + GUT.) ~ c^hc: ~ g 



-i(T,)G e -i(T e ) 
C+G e _j(Tj)+G 



aP + T (jujC e + G a + G e -i(T e ) 



G 2 



G e -i(Ti)G e -i(T e ) 

Ga+j^Ca juiC+Ge-tiTti+G 



.(19) 



2.3 Noise 



The effect of noise on a detector's performance is quantified by the Noise 
Equivalent Power (NEP). It corresponds to the power that would be required 
as input to the detector in order to generate an output equal to the signal 
generated by the noise. The NEP can therefore be calculated as the ratio 
between the output generated by the noise and the responsivity of the detector: 

A Y 

NEP, = ^ ,20) 



The variable y stands in for any of the possible noise terms: amj9=amplifier 
noise, j— Johnson noise, i?L=load resistor noise, 1//=1// noise, a=absorber- 
electron system thermal noise, th=h.e&t sink-lattice thermal noise or /ie=electron 
system-lattice thermal noise. 

To obtain the Noise Equivalent Power for each term, the noise contributions 
e amp , ej, e RL , P Rl , (AR/R) 1/f , P a , P th , and P he must be added to the block 
diagram of Fig. 3. Figure 4 shows where each noise term should be added to 
the block diagram (for more details see [3]). Solving the noise block diagram 
for each noise term independently and dividing by the responsivity obtained 
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in Eq. 16 we obtain the following NEP's: 

NEP a = P a (cu)jLUT a (21) 

NEP Rl = P RL {u){l+juT a ) + ^ (22) 

NEP amp = (23) 

NEP h = P h + (24) 

he ne (G + G e _,(T0+J^) (24) 

NEPth ~ ^(G + GUTd+juG) (25) 



NEP ej = ej (u)jj^(l +jur a )(G a + G e -i(T e )+juC e 

G e -i(Ti)G e -i(T e ) Gl 



juC + G^TD+G G a +juC a 



(26) 



NEP 1/f = +juT a )(G a + G e -i(T e ) + jujC e - 

Ge-ijTfiGe-ijTe) G\ \ 

juC + Ge-W+G G a +JiuC a )- 

Where r a = CJG a . 



3 Model 2 



In experiments involving dark matter detectors and double-beta decay detec- 
tors the absorber size is significant and can have a mass up to almost 1 Kg [5]. 
For mechanical reasons these large absorbers must be mechanically connected 
to the heat sink. The thermal link between the detector and the heat sink can 
therefore also be through the absorber rather than through the sensor. Model 
2 reflects this condition by having the absorber connected to the heat sink 
through a thermal conductivity G. Model 2 also takes into account absorber 
decoupling and electron decoupling by connecting the lattice system to the 
absorber through a thermal conductivity G a and by having the electron sys- 
tem connected to the lattice system through a thermal conductivity G e -i (see 
Fig. 5). 
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Applying the same procedure previously used for Model 1 we obtain for Model 
2 the block diagram of Fig. 6. Solving the block diagram for the responsivity, 
the dynamic impedance, and all the noise contributions, we obtain the follow- 
ing results: 



S(w)= 





Ga G e _i (Tjf ) 


XaAt r 








(G e -l(T e ) + G ETF +ju>C e ) 


(Ga+G^iT^+jwC^Ga+G+jwC^-Gl 


—G e _i (Ti)G e _i{Te)(G a +G+ju>Ca) 



(28) 



Z(u) = R- 



aP + T e 


jujC e + G e - 


l(Te) ~ 


-G e _ 


i(Ti){ 


G e _ i (T e )(jo;C a +G +G) \ 




(jujCa+Ga+G^jujd+Ga+Ge-lim-Gl J 




-aP + T e 


juiC e + G e 


-l(Te) 


-G e 


-i(Ti) 




)] 


\(3uCa+Ga+G){juC l +Ga+Ge-l{T l ))-Gl / 





(29) 



NEP th = P th 
NEP - 6amp 

1 ' ^ 1 amp 



NEP a = P a 



NEP he = P he 
NEP Rl = P Rl 



G + juiC a 



G a 



(30) 
(31) 

(32) 



(G a + G + juC a )(G a + G e -i(Ti) + juG) - G\ G + G a + juC a 



G a G e -i(Ti) 

(G a + G + juC a )(G a + G e -i(T t ) + juQ) - G\- 



G a G e -i(Ti) 



G n 



f ^(34) 



(33) 



S(u) 



NEP e = e j 

3 IRa 



Te (G a + G + juC a ) [(G a + G e _,(r,) + juQ) - l^gfejgf 1 ] - G\ 



G a G e -i(Ti) 



x(juC e + Gm) 



, AR(«>) \ (Ga + G + jujC a )[(G a + G e -i(ri) + juQ) - ^cgfejg 1 ] ~ ^ 

7 V Ji 



G a G e -i(Ti) 



x^(juC e + G e m 
a 



4 Examples of Energy Resolution and Time constant Results from 
the Models 



In experimental setups of detectors some of the parameters are fixed while 
others can vary. With the freedom to vary a few parameters the goal is to 
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optimize the detectors characteristics. The equations derived in this paper of- 
fer the flexibility and power to perform such operations. An example of such 
applications is reported in fig. 7. The characteristics of interest are energy- 
resolution and the time constant of the detector. The energy resolution is cal- 
culated following equation 69 in reference [3] . The time constant is calculated 
as the inverse of the first turn frequency of the responsivity. The variable pa- 
rameters in the examples of fig. 7 are: heat capacity of the absorber (C a ) vs. 
thermal conductivity between absorber and thermometer (G a ) and thermal 
conductivity between detector and heat sink (G) vs. thermal conductivity be- 
tween absorber and thermometer (G a ). The fixed parameters are reported in 
Table 1. By inputing the fixed parameter into the two models of this paper and 
the model found in [3] and letting C a , G , and G a vary we obtain the twelve 
contour plots in fig. 7. The first six plots predict energy resolution while the 
last six plots refer to the detectors time constant. Each column of contour 
plots belongs to one of the three different models. 



5 Conclusions 

To improve the performance of microcalorimeters and bolometers it is impor- 
tant to accurately understand how this depends on the fabrications param- 
eters. Significant improvements in detectors performance have, in fact, been 
achieved by optimizing the design based on an accurate model of the detector 
[6]. In this paper we derived detailed theoretical models to describe the behav- 
ior of two different detector architectures. The use of block diagram algebra 
has allowed us to present the results in an analytical form that can be easily 
and immediately utilized by investigators to improve the design of their detec- 
tors. Furthermore, the contour plots in fig. 7 provide an example of how the 
equations derived in this paper can be utilized to predict the characteristics 
of detectors. 



9 



References 



[1] J.C.Mather, Appl. Opt. 21, 1125 (1982). 

[2] S.H.Moseley, J.C. Mather, and D.McCammon, J. Appl. Phys 56, 1257 (1984). 

[3] M. Galeazzi and D. McCammon J. Appl. Phys. 93, 4856 (2003). 

[4] C. A. Kilbourne Nucl. Instrum. Methods. Phys. Res. A 520, 402 (2004). 

[5] C. Brofferio Nucl. Instrum. Methods. Phys. Res. A 520, 151 (2004). 

[6] M. Galeazzi, et al., Nucl. Instr. Meth. Phys. Res. A, 520, 469 (2004). 



10 



Tables 



Table 1: Values of the fixed parameters needed to produce the contour plots 
of fig. 7 



Figures 



FIG. 


1: 


Thermal architecture of Model 1. 


FIG. 2: 

biased, 


Typical readout circuit. Notice that if Rl « R the detector is voltage 
if Rl » R the detector is current biased. 


FIG. 


3: 


Block diagram representing Model 1. 


FIG. 


4: 


Block diagram including noise contributions for Model 1. 


FIG. 


5: 


Thermal architecture of Model 2. 


FIG. 


6: 


Block diagram including noise contributions for Model 2. 


FIG. 


7: 


Contour plots of how the energy resolution and time constant of each 



model change with respect to the heat capacity C a and the thermal conduc- 
tivities G and G a . These plots were constructed using the fixed parameters in 
Table 1. 
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Table 1: Values of the fixed parameters needed to produce the con- 
tour plots in fig 7 



Parameter 


Value 


R 


5mQ 


Rl 


0.2mQ 


Vbias 


1.0848 x 10- 7 V 


a 


100 


c a 


IpJ/K 


Q 


4.911 x 10~ 5 pJ/K 


c e 


0.15A P J/K 


T s 


0.1K 




5 x 10~ W W/K 


G e -l(T e ) 


5.87 x 10- 10 W/K 


G 


1 x 1Q- 10 W/K 



12 



w 

C a v/sAJ C e i 



G e _i 

c, 

> G 
'/////>////. 

HEAT SINK 

Fig. 1. Thermal architecture of Model 1. 




Fig. 2. Typical readout circuit. Notice that if Rl « R the detector is voltage 
biased, if Rl >> R the detector is current biased. 



AT, (o>) 



G + G_ 1 (T,) + jcoC 1 — G e-i(T e ) ~| 



W(co) 



1 1 



lAT a (ffl) 



G a l+jcox a 




G e4 (T e ) + G a +jcoC e 



AT e (<D) 



XaA 



Gf 



AX(co) 



G s 



Fig. 3. Block diagram representing Model 1. 
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p th («» 



AT, (a)) 



W(a>) 

-ty" 



1 



G a 1+jfflT, 



AT a (co) 



G + G„ 



ej(») 



P he^ 

P Rl V) 



1 



Ge-l( T e) + G a +j»Ce 



89 



T 

IK a 



&6 



AT e ((B) 



X-a-A, 



e a mp( ffi ) 

+ 

eR L («>) 
■J*AX(o>) 



Fig. 4. Block diagram including noise contributions for Model 1. 



P=I 2 R 



W 



*e-l 



G a 



777777777, 

HEAT SINK 



Fig. 5. Thermal architecture of Model 2. 



P th (a» 



w(a>) I i 



'G. + G 1 + jffl 



-Tg, 



^ G e-lt T e)r 



G„ I c 



G,+G e ,(T,) 



1 + jO) 



Ga+GelCT,) 



ATi(o) 



G e-l( T eH + j G)T e-) 



AT,.(m) 



(AR1 

lRJ 1/f < m > «amp(») 



E« a I 



Fig. 6. Block diagram including noise contributions for Model 2. 
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Galeazzi & McCammon 2003 Model 1 (This Paper) 

Energy Resolution Go vs Co Energy Resolution Go vs Co 



Co (W/K) 



Energy Resolution Ga vs G 




Time Constant Ga vs G 




Model 2 (This Paper) 

Energy Resolution Ga vs Co 



Ga (W/K) 

Energy Resolution Ga vs G 



Go (W/K) 
Time Constant Ga vs Ca 



Ga (W/K) 
Time Constant Ga vs G 




Time Constant Ga vs G 



10 -id 

Go (W/K) 




Go (W/K) 




0.1 1.0 10.0 100.0 0.0001 0.0010 0.0100 0.1000 

Energy Resolution (eV) Time Constant (s) 



Fig. 7. Contour plots of how the energy resolution and time constant of each model 
change with respect to the heat capacity C a and the thermal conductivities G and 
G a . These plots were constructed using the fixed parameters in Table 1. 
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